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SPECIFICATION 

[Electronic Version 1.2.8] 

METHOD AND SYSTEM FOR TIME 
DIFFERENCE OF ARRIVAL (TDOA) 
LOCATION SERVICES 

Background of Invention 

[0001] Field of the Invention. This invention relates to methods and systems for 

location finding. More specifically this invention relates to location finding methods and 
systems, which make use of "Time Difference of Arrival." Description of Related Art. A 
variety of methods and techniques have been proposed for geographic location 
determination. The most widely used system is the Global Positioning System (GPS). 
The GPS system is a satellite-based pseudo-ranging location system that is used widely 
in military and commercial applications. Typical GPS accuracy is less than 10 meters for 
military applications and 100 meters for commercial applications. Differential GPS 
(DGPS) has been introduced to improve the location accuracy for commercial 
applications. GPS provides accurate position information for users in open and 
unobstructed areas. It requires that at least three satellites be within range of the GPS 
receiver to provide 2-D position data and four or more satellites to provide 3-D data. 
Satellite blockage caused by man made or natural structures within urban and/or 
shadowed environments can render GPS inoperable. Generally such GPS systems do not 
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provide the broad outdoor coverage in an urban and/or shadowed environment. Also, 
such GPS systems tend to be relatively costly. 

. [0002] Other common position location methods include Direction Finding (DF), 

hyperbolic methods and Time of Arrival (TDOA). DF techniques generally require 
steerable antennas or phased-array antenna systems. Although DF techniques can offer 
substantial spatial resolution and have the ability to locate a number of signals, typically 
DF systems are highly complex. The complexity translates to high cost, low reliability 
and often imposes size and power restrictions on the use of the system. While TOA is a 
simple triangulation method, it requires both the base stations and all of the mobile 
devices to have highly accurate clocks. This requirement for highly accurate clocks 
imposes a significant cost impact on the use of the TOA method. 
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[0003] Although the following may not qualify as "prior art," the reader is referred to 

the following U.S. patent documents for general background material. Each of these 
referenced patents is hereby incorporated by reference in its entirety for the material 
contained therein. U.S. Patent Nos.: 3,996,590; 4,060,809; 4,252,999; 4,783,744; 
5,689,269; 5,752,218; 5,781,150; 5,952,959; 5,963,130; 6,041,222; 6,061,337; 
6,104,978; 6,184,801; 6,198,390; 6,204,81 3; 6,236,365; 6,252,544; 6,249,252; 
6,295,023; 6,324,21 3; 6,353,406; 6,388,61 1 ; 6,392,692; 6,424,826; 6,429,81 1 ; 
6,429,81 2; 6,452,544; 6,459,704; 6,497,656; 6,509,829; and 6,560,462. 
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Summary of Invention 



[0004] It is desirable to provide a method and system for locating a telemetry 

communication source position, while minimizing location error. More specifically, it is 
desirable to provide a method and system for identifying three dimensional location 
coordinates, using a Time Difference of Arrival (TDOA) method. 

[0005] Accordingly, it is an object of this invention to provide a method and system 

for determining a geographic coordinate location that is appropriate for use in a 
telemetry communication system. 

■' * . 

[0006] Another object of this invention is to provide a method and system for 

determining a geographic coordinate location that provides broad outdoor coverage in a 
shadowed environment. 

[0007] A further object of this invention is to provide a method and system for 

determining a geographic coordinate location that does not require the clocks at the 
receiver sites of interest to be synchronized to the transmitter clock. 

[0008] A still further object of this invention is to provide a method and system for 

determining a geographic coordinate location that employs a triangulation technique 
based on Time Distance of Arrival (TDOA). 
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[0009] Additional objects, advantages, and other novel features of this invention will 

be set forth in part in the description that follows and in part will become apparent to 
those of ordinary skill in the art upon examination of the following, or may be learned 
with the practice of the invention as described herein. The objects and advantages of 
this invention may realized and attained by means of the instrumentalities and 
combinations particularly pointed out in the appended claims. Still other objects of the 
present invention will become readily apparent to those skilled in the art from the 
following description wherein there is shown and described the preferred embodiment 
of the invention, simply by way of illustration of one of the modes best suited to carry 
out this invention. As it will be realized, this invention is capable of other different 
embodiments, and its several details and specific steps are capable of modification in 
various aspects without departing from the invention. Accordingly, the objects, 
drawings and descriptions should be regarded as illustrative in nature and not as 
restrictive. 

Brief Description of Drawings 

[001 0] The accompanying drawings incorporated in and forming a part of the 

specification, illustrate present preferred embodiments of the invention. Some, 
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although not all alternative embodiments are described in the following description. In 
the drawings: 

[001 1] Figure 1 is a system block diagram showing the major components of one 

preferred embodiment of the invention. 

[001 2] Figure 2 is a process flow chart of one preferred embodiment of the method 

of this invention. 

[001 3] Reference will now be made in detail to the present embodiments of the 

invention, examples of which are illustrated in the accompanying drawings. 

Detailed Description 

[001 4] This invention is a method and system for the locating of antenna device that 

uses a Time Difference of Arrival (TDOA) method. TDOA is a technique that does not 
require the mobile device that is to be located, to have an incorporated clock. Rather 
. than having known Euclidian distances of the Time of Arrival technique, TDOA uses the 
range differences between receivers. These range differences can be described as a 
hyperbolic curve in 2D or a hyperboloid in 3D space. This hyperbolic position location 
system estimates the location of a source by the intersection of hyperboloids describing 
range difference measurements between three or more base stations. For example; in 
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the case of a 2D hyperbolic solution to the location finding problem with three base 

station sites (SI through S3) and a mobile transmitting device (source) to be located. 

The relationship between the range differences and the TDOA between the receivers is 

given by: 
* 

[0015] Rij = cxdij = Ri-Rj 

[001 6] where c is the signal propagation speed (speed of light = 3 x 1 0 8 

meters/second), d(ij) is the TDOA between receiver i and j. The TDOA estimate, in the 
absence of noise and interference, restricts the possible source locations to a 
hyperboloid of revolution with the receiver as the foci. 

[001 7] In a 3D system, the hyperboloids that describe the range differences, R(i j) 

between receivers are given by: 

[001 8] Ru = [(Xi -x)2 + (Yi - y)2 + (Zi - z)*yn - [( Xj -x)2 + (Yj - y)2 + (Zj - z)2]i/2 

[001 9] where [(X(i), Y(i), Z(i)] and [X(j), Ytf), Z(j)] define the location of the receiver i 

and j respectively. R(i j) is the range difference measurement between base station i and 
j, and (x,y,z) are the unknown source coordinates. 

[0020] For the normal solution case, if the number of unknowns or coordinates of 

the source to' be located is equal to the number of equations or range difference 

» 
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measurements, then the system is consistent and a unique solution exists. Therefore, 
for a 2D system, where the source and the receivers are coplanar, a two-dimensional 
source location is determined from the intersection of two or more hyperbolas produced 
from three or more TDOA measurements. While, for a 3D system, the intersection of 
three or more independently produced hyperboloids generated from four or more TDOA 
measurements are used. 

[002 1 ] The triangulation system of this invention uses this TDOA technique. In the 

present embodiment of this invention, each base station site is provided with a Global 
Positioning System (GPS) receiver that is used to synchronize the towers to the cesium 
clock standard used in the GPS system and also provides a 1 0 nano-second (1 00 MHz) 
reference clock that is used in the measurement of the TDOA. Alternative embodiments 
may uses different and improved GPS systems and can employ different reference clocks 
without departing from the concept of this invention. This present approach has an 
inherent system error of plus or minus 1 0 nsec, which is approximately 1 foot of 
distance error per nano-second, or about 10 feet of location determination error. 
Multipath error is the dominant error source in this system. 

[0022] The solution to the TDOA location finding of this invention uses essentially 

the same set of equations as in the above described spherical case, but typically the R 
values are also unknowns. As an example, the following MathCad model, which solves 
the case for three receiver sites and one particular device position, is shown: 
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[0023] mi:= 5280; set 1 mile equal to 5280 feet 

[0024] site 1 := (0, 0, 0.1) x mi; set positions of towers [1 , 2, 3] 

[0025] site 2:= (5, 0, 0.2) x mi 

[0026] site 3:= (5, 5, 0.3) x mi • 

[0027] noise:= 10; position of device in feet 

[0028] dev:= (3, 4, 1) x mi 

t 

[0029] dev:= (1 .584 x 1 0<, 2.11 2 x 1 0\ 5.28 x 1 0*) 

[0030] Euclid(a,b):= [(a 0 - b 0 ) 2 + (ai - bi) 2 +' (a 2 - b 2 ) 2 ] 1/2 ; generate function to 
return 

[0031] ; Euclidian distance between 2 

[0032] ; xyz vectors . 

[0033] c:= (1860005280) c= 9.821 x 10*; speed of light in feet / second 
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[0034] E:= [Euclid(sitel ,dev) + rnd(noise), ; reported position from towers and 
noise 

[0035] Euclid(site2,dev) + rnd(noise) f Euclid(site3,dev) + rnd(noise)] 

[0036] E:=[2.682 10«, 2. 399 xlOM .238 x 10*] 

[0037] T:= 100; set realtime stamp to 100 

[0038] T:= T + E; simply so RF travels in unit time, therefore 

[0039] T = (2.692 x 1 0", 2.409 x 1 0", 1 .248 x 1 Of) ; time at each site is. T + E, actual 

[0040] ; formula is T = T + E/c, but the huge 

[0041] tl = 2.409 x 10 4 ; c value causes scaling problems for 

[0042] ; the solver 

[0043] delta 0 := (T 0 - T,); We also can calculate the delta times between 

[0044] deltai:= (To - T 2 ); the stations 
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[0045] delta 2 := (T, - T 2 ) 



[0046] delta = (2.835 x 103, 1.445 x 10\ 1.161 x 10<) 



[0047] u:= (3, 3 ? .5) x mile R:= (3, 3, 3) x mile; set guess values for solver 



[0048] Given; device report z with barometer, 



[0049] ; therefore u[z] is a constraint 



[0050] R 0 - Ri = deltao 



[0051] Ro-R2 = delta, 

Ri - R2 = delta2 u 2 = dev 2 ; set z to baro altitude 



[0052] A system of simultaneous equations, given known xyz's of sites, verses 

unknown (u) xyz's of device to be located, using received Euclidean distances R are 
solved. 



[0053] (sitel 0 - u 0 ) 2 + (sitel 1 - uO 2 + (sitel 2 - u 2 ) 2 = Ro 2 



[0054] (site2 0 - u 0 ) 2 + (site2, - u0 2 + (site2 2 - u 2 ) 2 = Ri 2 



[0055] (site3 0 - u 0 ) 2 + (site3i - U1) 2 + (site3 2 - u 2 ) 2 = R2 2 
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[0056] Triangulate(u,R):= Minerr(u.R) 

[0057] ans:= Trangulate(u.R) 

[0058] cdev:= ansO 

[0059] cdev = (1.584 x 10", 2.112.x 10", 5.28 x lO*) 

[0060] dev = (1.584 x 104,2.112 x 10", 5.28 x 103) ' 

[0061] cdev- dev = (-2.075, -2.34, 0) 

[0062] Figure 1 shows a system block diagram showing the major components of 



one preferred embodiment of the invention. A satellite system 100 is used to facilitate 
and the operation of the GPS receiver 102a-d located in the base stations lOla-d. the 
GPS receiver 102a-d locks the base. station lOla-d clocks 105a-d to the NTSB Cesium 
clock standard. In the present embodiment of the invention, both the carrier frequency 
as well as the A/D sampling clocks of the base station lOla-d are locked and derived 
from a 10 MHz signal provided by the GPS receiver. Each base station lOla-d is 
provided with a GPS receiver 102a-d, a processing unit 103a-d, an altitude (typically 
barometric) sensor 104a-d and, in the present embodiment, a reference clock 105a-d. 
The GPS receivers 1 02a-d are preferably provided with an antenna 106a-d, that is in 
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electronic communication with the GPS receiver 1 02a-d. The CPS receiver 1 02a-d, the 
altitude sensor 104a-d and the reference clock 105a-d are typically in electric 
communication with and, to varying degrees, controlled by the processor 1 03a-d. 

[0063] Figure 2 shows a process flow chart of one preferred embodiment of the 

method of this invention. If the signal received from the device to be located is received 
by four or more receiving sites, the TDOA triangulation method as follows is performed. 
Alternatively, if the signal received from the device to be located is received only by 
three receiving sites, an assumption is initially made that the three sites and the device 
are all in the same plane. The z-value of the solution to the mean fixed altitude of all* 
three sites is constrained and the following TDOA triangulation method is performed as 
follows. This approach obtains an approximate position of the device to be located., 
Then a surface interpolation of each site's barometric altitude is performed and the 
difference in altitude between the device's reported barometric altitude and the site's 
interpolated barometric altitude is used to compute an actual z-altitude of the device to 
be located, based oh the known z-altitudes of the sites. The altitude estimate of the 
device to be located is used to constrain the z-value of the site (which did not receive a 
signal from the device) having the unknown xyz position and the TDOA triangulation 
method, as follows, is performed. This approach provides a resulting xyz position that 
is the current best estimate of the position of the device to be located. 
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[0064] The present TDOA triangulation process of this invention is detailed in figure 

2. The initial values are set 201 . These initial values typically include such constants as 
feet to miles and the positions of the towers (base station sites). The Euclidian 
distances between the site vectors are calculated 202. The reported position is 
calculated 203 with noise. The realtime stamp is set 204. The time stamp is simplified 
204 so that scaling problems, not critical to location determination, are removed. The 
delta times are calculated 206 between the stations. A guess value for the solver is set 
207. The z value is reported 208 from the altitude sensor (typically a barometer). The 
Euclidean distances of the sites and unknown devices are calculated 209. The location 
of the unknown device is then triangulated 210. This triangulation 210 step can include 
the generation of deviation statistics and error correction as further described as 
follows. 

[0065] Although the range differences calculated from the data provided by each 

site lOla-d can drift over the combined circular error probability (CEP) of the GPS 
receivers 102a-d. The resulting error in the range difference data is corrected, in the 
present embodiment of the invention, by (1) determining the actual xyz location of each 
site from map or survey data; (2) calculating the error between the reported xyz position 
of the base station lOla-d from the base station GPS receiver 102a-d; (3) collecting 
range difference data from xyz positions from each GPS receiver 1 02a-d relative to the 
known locations of the base station sites; and (4) use the range difference data, in the 
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simultaneous equations with the base station site locations set from the known xyz 
locations to converge on the correct location position. 

[0066] The described embodiment of this invention is to be considered in all 

respects only as illustrative and not as restrictive. Although specific steps and 
associated formulas are provided, the invention is not limited thereto. The scope of this 
invention is, therefore, indicated by the claims rather than by the foregoing description. 
All changes, which come within the meaning and range of equivalency of the claims, are 
to be embraced within their scope. 
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